Constraining the variation of fundamental constants using 18cm OH lines 
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We describe a new technique to estimate variations in the fundamental constants using 18cm OH 
absorption lines. This has the advantage that all lines arise in the same species, allowing a clean 
comparison between the measured redshifts. In conjunction with one additional transition (for 
example, an HCO + line), it is possible to simultaneously measure changes in a, g v and y = m e /m p . 
At present, only the 1665 MHz and 1667 MHz lines have been detected at cosmological distances; we 
use these line redshifts in conjunction with those of HI 21cm and mm-wave molecular absorption in 
a gravitational lens at z ~ 0.68 to constrain changes in the above three parameters over the redshift 
range < z < 0.68. While the constraints are relatively weak (< 1 part in 10 3 ), this is the first 
simultaneous constraint on the variation of all three parameters. We also demonstrate that either 
one (or more) of a, g p and y must vary with cosmological time or there must be systematic velocity 
offsets between the OH, HCO + and HI absorbing clouds. 

PACS numbers: 98.80.-k,98.80.Es,98.58.-w,33.15.Pw 



INTRODUCTION 

The recent claim by Webb et al. 0,0 that the fine structure constant a evolves with redshift, with Aa/a = 
(-1.88 ± 0.53) x lCT 5 from z ~ 1.6 to today and Aa/a = (-0.72 ± 0.18) x 10~ 5 for 0.5 < z < 3.5 (but see ;3]) has 
spurred interest in the possibility that the numerical values of the fundamental constants change with time. Theories 
that can account for such variations include extra-dimensional Kaluza-Klein theories and superstring theories. In 
such models, the values of the coupling constants depend on the expectation values of some cosmological scalar 
field(s); changes in the values of the coupling constants are thus to be expected if this field varies with location and 
time. Further, depending on the details of the theory, all of the different coupling constants (such as a, the proton 
g-factor g p , the electron-proton mass ratio y = mJm v , the gravitational constant G, etc) could, in principle, vary 
simultaneously. For example, Calmet & Fritzsch £| and Langacker et al. 5] find that variations in the value of a 
should be accompanied by much larger changes (by ~ 2 orders of magnitude) in the value of y. However, Ivanchik et 
al. [|| constrain the variation in y to be (3.0 ±2.4) x 10~ 5 over the redshift range < z < 3, comparable to the change 
claimed in the fine structure constant. A review of the available experimental and observational measurements on the 
variability of the coupling constants can be found in 

One of the main problems in most of the astrophysical techniques used to measure (or constrain) the values of the 
different constants (e.g. @,E3,@) i s that they involve a comparison between the redshifts of spectral lines of different 
species (e.g. the HI 21cm line, millimetre-wave molecular lines and optical fine structure lines Hi). These species 
are unlikely to all arise at the same physical location in a gas cloud and might thus have systematic velocity offsets 
relative to each other; the redshift differences may thus be dominated by these effects rather than the measurement 
errors (i.e. the spectral resolution, which can be quite small, ~ few km/s, for HI 21cm and mm-wave molecular 
absorption spectra). Conclusions drawn from a comparison between different species might thus well be incorrect. 

Clearly the best way to test the variation of the coupling constants is to use lines originating from a single species, 
but with different dependences on these constants. Further, since many constants may be varying simultaneously, 
it would be very useful if one could simultaneously measure the changes in a number of constants from this single 
species, rather than assuming that changes occur in only one of the constants and that the others remain unchanged. 
We present here a new technique which satisfies both these requirements, using the 18cm lines of the OH radical. 

The ground 2 n 3 / 2 J=3/2 state of OH is split into two levels by A-doubling and each of these A-doubled levels further 
split into two hyperfine states. Transitions between these levels lead to four spectral lines with wavelength ~ 18cm. 
Transitions with AF = are called the "main" lines, arising at rest frequencies of 1665.4018 MHz and 1667.3590 MHz, 
while transitions with AF = 1 are called "satellite" lines, with rest frequencies of 1612.2310 MHz and 1720.5299 MHz. 
Since the four OH lines arise from two very different physical processes, viz. A-doubling and hyperfine splitting, the 
transition frequencies have different dependences on the fundamental constants. A perturbative treatment of the OH 
molecule has been carried out by Dousmanis et al. [!j (see also ^(J); we use the expressions derived in these references 
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to determine the dependence of linear combinations of the four OH line frequencies on a, g p and y = m e /m p and show 
that it is possible to simultaneously measure changes in both a and y, if all four line frequencies are known (assuming 
that g p does not vary with time). Since OH and HCO + column densities are observed to be tightly correlated, both 
in the Galaxy and out to z ~ 1 [l!^. these species are likely to arise at the same physical location; a comparison 
between the redshifts of the four 18cm OH lines and the HCO + line should thus allow one to constrain the evolution 
of all the above three parameters. We use our observations of the OH main lines in the z ~ 0.6846 gravitational 
lens towards B0218+357, in tandem with published HI 21cm and mm-wave molecular redshifts, to constrain Ay/y 
between z ~ 0.68 and today. Finally, as this work was being written up, an analysis on the use of OH lines to constrain 
changes in fundamental constants was also carried out by Darling [lfgj ; the latter, however, only considers variations 
in the fine structure constant a. 



CONSTRAINTS FROM RADIO SPECTRAL LINES 



Consider two transitions whose rest frame frequencies vi(z\i — 1,2 depend on redshift, due to the evolution of 
various fundamental constants, such as a, g p ,m e /m p , etc. If the lines arise in a source at a "true" redshift z, the 
measured redshift Zi of each line is given by 
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where Avi(z) = Vi(z) — Vi(0); note that Avi(z) = if the fundamental constants do not change with time. The first 
order difference between the two measured redshifts Az = z\ — z^ is then 
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where z is the mean measured redshift. Given two spectral lines (or linear combinations of line frequencies) with 
different dependences on some fundamental parameter, one can thus use the differences between the measured redshifts 
to constrain the evolution of the parameter in question. 

In the case of the four 18cm OH lines, the following three independent relations have been shown to be satisfied by 
the line frequencies |jl Il0|: note that the lines must also satisfy the constraint i^i665 + ^1667 = ^mo + ^1612 
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Equations © , & JSJ correspond to the energy split due to A doubling and the difference and sum of the hyperfine 

1/2 

splits in the two A doubled levels respectively. Here, X = [(A/B) {(A/B) — 4} + 16] , A is the fine structure 
interaction constant, B, the rotational constant, A' and B' , the off-diagonal matrix elements of these operators, q\ ~ 
4B 2 /his e (where hv e is the energy difference between the ground and first excited electronic state), and, finally, a, 6, c 
and d are "hyperfine constants" 9], whose experimental values are a = 86.012±0.002 MHz, b = — 116.719±0.008 MHz, 
c = 130.75 ± 0.01 MHz and d = 56.632 ± 0.004 MHz 14]. Numerically, A/B = -7.547 and A'/B' = -6.073 [l^. 
These quantities have the following dependences on the fundamental constants a, y, Roo = m e e^/h 3 c and g p : 
A 1 cx A oc a 2 Rcyr, , B' oc B oc yRoo and a, 6, c, d oc gpO^yRoc [lOj. For the rotational constant B, we have assumed, 
as usual (e.g. [l5|)) that variations in (m p /M), which are suppressed by a factor m p /U ~ 100 (where M is the reduced 
mass and U the binding energy) can be ignored. Thus, we have [A'/B'] oc [^4/B] oc (a 2 /y). Replacing the above 
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scalings in equation @ for z/4, we obtain z/4 oc y 2 RooF {a 2 /y) , where F = F (/3) is a function which depends only 
on the ratio (3 = A/B oca 2 jy and is defined by 
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Evaluating the quantity on the right-hand-side of the above equation, we obtain 

a Ay Aa AR^ 
= 2.571— - 1.141 + - 

v A y a Rao 



In similar fashion, equations Q and JSJ yield 
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Equations (0 - 111(1 govern the way in which a change in one of the fundamental constants affects the line rest 
frequencies. Combining them in pairs in equation @ yields : 
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If all four OH lines are detected in absorption in a single cosmological system, we thus have two independent equations 
relating the differences in measured redshifts to the changes in a, g p and y = (m P lm n ). Most analyses (using other 
lines) assume that g p remains unchanged and then estimate the variation of a 0, llof . If one makes the same 
assumption in the case of the OH lines, equations (|12I - I13(I immediately allow one to simultaneously solve for changes 
in both a and (m e /m p ). Of course, one further equation is needed to simultaneously constrain the evolution of all 
three constants. Candidates include the HI 21cm and mm-wave molecular lines. The best of these are likely to be 
the HCO + lines since HCO + and OH column densities are found to show a strong correlation (extending over more 
than two orders of magnitude in column density) both in the Galaxy ^l| and out to z ~ 1; this suggests that 
HCO + and OH are located in the same region of the molecular cloud. Since the HCO + line arises from a rotational 
transition, we have 
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Equations (|12fl . I|13l) and 1141) all have the same dependence on the Rydberg constant Roo, which hence cancels out 
(note that Roo itself depends on a through the relation Roo = m e e 4 /h 3 c). One might thus combine HCO + absorption 
redshifts with those derived from the 18cm OH lines to provide the last equation needed to solve for the evolution of 
a, y and g p . This would allow a simultaneous measurement of all three constants, which has been hitherto impossible. 
We note that it is also possible to use other OH A-doubled transitions to simultaneously constrain the variation of 
the above parameters; this is discussed elsewhere ^||, for the "main" OH lines. It would also be interesting to carry 
out an analysis similar to that of Bekenstein Q to test whether these analyses of the OH lines are also affected by the 
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possibility that the Hamiltonians involved might vary with time (i.e. dynamical variability of the different parameters); 
this is beyond the scope of the present paper, ft should also be pointed out that the present calculation is based on 
an analysis of the OH levels using perturbation theory Q; more recent analyses 01 use the "effective Hamiltonian" 
approach, giving rise to higher order effects. While it would be interesting to attempt the present calculation in the 
latter framework, we emphasise that it should only result in small changes to the coefficients in equations (Pl- lllf) 
and does not affect the validity of our approach. Finally, while other "Lambda-doubled" systems could, in principle, 
be used for a similar analysis, none of these have multiple transitions detected in astrophysical objects (to the best of 
our knowledge). While it would be interesting to carry out searches for these other transitions, we suspect that OH 
is yet likely to prove the best candidate because of the strength of its multiple lines. 



APPLICATION TO THE z = 0.6846 ABSORBER TOWARDS B0218+357 

While the above analysis shows that one can use the four f 8cm OH lines (in conjunction with the HCO + transition) 
to constrain the evolution of three separate fundamental parameters, the weaker, satellite 1612 MHz and 1720 MHz 
lines have so far not been detected at cosmological distances. Observations are currently being scheduled to carry out 
deep searches for these lines in the four known OH absorbers at intermediate redshift [lil [Lsl ITfl) . For the present, 
we will instead use the detected 1665 MHz and 1667 MHz transitions (i.e. equations © and l|lQ|l ) along with the 
HI 21cm and millimetre (HCO + ) lines in the z = 0.6846 absorber towards B0218+357 to estimate changes in the 
fundamental parameters (assuming that the lines arise in the same gas cloud). Since the HI 21cm frequency arises 
from a hyperfme split and is hence proportional to g P ycx 2 R 00 , we have 

Ais 21 _Ay | /a | Ag p | AR^ 
vi\ y a g p 

Equations Q, l|10(l . (|14Jl and 1)15(1 can now be solved to measure changes in a, g p and m e /m p . The HI 21cm redshift is 
Zhi = 0.684676 ±0.000005 ||, while that of the HCO+ absorption lines is z HCO + = 0.684693 ±0.000001 20]. Our new 
GMRT OH absorption spectra towards B0218±357 yield the following redshifts for the sum and difference of the 
1665 MHz and 1667 MHz line frequencies : z sum = 0.684682±0.0000056 and z diff = 0.685780±0.0067. A simultaneous 
solution of equations ©, O, 03 and EH then yields (Aa/a) = (-0.38±2.2) x 10~ 3 , (Ay/y) = (-0.27±1.6) x 10" 3 
and (Ag p /g p ) = (—0.77 ±4.2) x 10~ 3 . Since the error on Zfctt is far higher than the other errors, this dominates the 
errors on the above estimates and results in relatively uninteresting upper limits on changes in the three constants. 
We emphasise, however, that, to the best of our knowledge, this is the first time that a simultaneous constraint on 
the variation of these three fundamental parameters has been obtained in a cosmologically distant object. Further, 
the weakness of the constraint arises from the relatively small difference between the frequencies of the main lines; 
this is clearly not a fundamental limitation but depends entirely on the sensitivity of the observations. 

Next, if we assume that y = m e /m p is constant, equations © and ITljl yield [Aa/a] = —5.7 ± 3.0 x 10~ 6 . This is 
more than 3er deviant from the estimate [2 Aa/a] = 1 ±0.3 x 10~ 5 of Carilli et al. (2000) ||. Since the latter analysis 
assumed that g p was constant, the difference between the two estimates implies either that the assumptions that y 
and/or g p are constant is unjustified or that systematic velocity offsets do exist between the three species. 

Further progress can be made if one of the three parameters, a, g p and y is assumed to not change with time 
(while retaining the assumption that velocity offsets are not significant). We can then avoid having to use the 
equation governing the difference between the OH redshifts and can thus obtain a far stronger limit on changes in 
the remaining two quantities. For example, if we assume (as is often done, e.g. |8») that g p remains constant, the 
HI 21cm and mm-wave molecular line redshifts imply [Aa/a] = (5 ± 1.5) x 10 -6 Combining equations © and 
(I14JI then yields [Ay/y] — (7.8 ± 2.4) x 10~ 6 . Similarly, if we assume that y = m e /m p is constant, equations © and 
<|Tif yield [Aa/a] = (-5.7 ± 3.0) x 10" 6 and [Ag p /g p ] = (2.2 ± 0.67) x 10~ 6 . Finally, if a remains unchanged, we 
obtain [Ay/y] = (4.2 ± 2.2) x 10~ 6 and [Ag p /g p ] = (1 ± 0.3) x 10" 5 . The above limits on Ay/y are a factor of - 4 
stronger than the best earlier limits on changes in this quantity [6j. It is very interesting that all three cases result 
in a higher than 3tr significance for the variation of at least one of the parameters. This implies either that one (or 
more) of these parameters indeed varies with cosmological time or that systematic motions between the three species 
cause the above uncertainties (which only include measurement errors) to be under-estimated. 

Finally, we note that, while the OH 1667 MHz and 1665 MHz, HCO + and HI lines have been detected in four 
absorbers at intermediate redshifts, two of the absorbers (PKS1413±135 and B2 1504±377) are believed to have 
velocity offsets between the HI and HCO+ redshifts [EEl. The last system, at z ~ 0.889 towards PKS1830-21, does 
not at present have OH data of sufficiently high quality to carry out the above analysis. 
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In summary, we have demonstrated a new technique to simultaneously measure the evolution of the three funda- 
mental constants a, g p and m e /m p , using 18cm OH absorption lines in conjunction with one additional transition, 
(which could be an HCO+ mm- wave line). At present, only the 1665 MHz and 1667 MHz "main" OH lines have been 
discovered at cosmological distances; we have used these line redshifts in conjunction with those of HI 21cm absorption 
and millimetre-wave molecular lines to constrain the variation of y, g p and a between z — 0.6846 and today. We argue 
that one (or more) of the parameters a, y and g p must vary with cosmological time, unless systematic velocity offsets 
exist between the above three species. The constraints placed on changes in the parameter y = m e /m p (assuming 
that cither a or g p are constant) are a factor of ~ 4 stronger than earlier limits on variations in this parameter. 

We are grateful to Rajaram Nityananda for very useful discussions on the energy levels of the OH ground state. 
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